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Abstract 

We calculate the energy of electroweak sphalerons including one-loop fermionic correc- 
tions. This calculation has previously been done by adding the correction to the tree-level 
sphaleron and interpreting the resulting energy as the 1-loop energy. However, in this paper 
we calculate the energy by doing a full renormalisation of the parameters on 1-loop-level 
and redetermining the sphaleron configuration using the full energy functional. When com- 
paring the final result with the tree-level solution we find that the 1-loop calculation will 
only cause small deviations of the sphaleron energy. 
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1 Introduction 



In the study of non-abelian gauge theories it has been found that the non-trivial 
structure of the gauge fields can lead to baryon number violating processes due to 
the anomaly of baryon and lepton currents (see |]J). This fascinating property of 
Yang-Mills theory was discovered by Faddeev and by Jackiw and Rebbi who 
showed that different field configurations of the gauge-fields exist that are topologi- 
cally distinct but are physically all vacuum configurations. Thus the field energy of 
these configurations will be zero. These configurations are numbered by a certain 
functional of the fields, the Chern-Simons number iVcs (which is, mathematically 
speaking, the winding number). They can be continuously deformed into each other 
but this process will violate the boundary condition that the energy is zero. However, 
there exists a way of deforming the fields that will make the required energy as low 
as possible. The peak of this path is known as the sphaleron and its energy is of vital 
interest for the whole process. The energy of the sphaleron is of the order of m^ja 
where mw is the mass of the W boson and a the SU{2) coupling constant. 
The classical sphaleron energy (the tree-level result) has been known for some time 
already. Naturally the question arises, how accurate the tree-level result is. Quantum 
corrections to the sphaleron have been computed by Bockharev and Shaposhnikov; 
they have included boson fluctuations through an effective potential of the Higgs field 
(see [f|, [5]]). A direct computation of the bosonic determinant over nonzero modes has 
been performed by Carson and McLerran || [7| and Baacke and Junker |§. 
In || for the first time fermionic quantum corrections both at zero and at finite 
temperature were examined and their influence turned out to be significant. However, 
these calculations performed the renormalisation only approximately and they did not 
take into account that the gauge fields might look different on 1-loop level than on tree- 
level. The generalisation of these result to the case including bosonic contributions 
was done in [[Kj]. However, the conclusions drawn there were based on the above- 
mentioned simplifications in the renormalization and differ from the results given by 
Moore in [[I I]] and Shaposhnikov [12j . 

The computation of the one-loop sphaleron configuration by using a self-consistent 
approach is the primary aim of our paper. We further perform a correct renormali- 
sation of all parameters of the theory and discuss the effects on the results given in 

H- 

In section 2 we discuss the background of sphaleron configurations and fermionic 
fluctuations. Section 3 discusses the methods and results of the self-consistent cal- 
culation. In the last chapter we draw the conclusions and relate this work to other 
papers. 
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2 Sphalerons and fermionic fluctuations 



The particle model under examination is the minimal standard model with one Higgs 
doublet. In |l3j it has been analysed, how large the influence of the {/(l)-sector to 



the sphaleron energy is. As its influence is only of the order of ~ 1% it seems to 
be justified to neglect the {/(l)-sector in the model (by setting the Weinberg angle 
$w = 0) an d work only with the S'?7(2)-symmetric Yang-Mills theory with a single 
Higgs doublet. The Lagrangian is thus (for one fermion species) 



v2 / „,2\ 2 



-^ F U Fa " v + (pM P"*) - y - yj (!) 

where we introduced the covariant derivative _D M = d^ — iA^ with = ^A^T a and the 
field strength tensor F^ = \F^ v r a , where F£ v = d^A a v - d u A° + e ahc A b ^A c v . The mass 

matrix M consists of the components of the Higgs field $ = {^o \ an d the Yukawa 
couplings h u and h d , 



o* 



The SU (2) fermion doublets are defined by 

R/ 



2 V ' / r 

Masses are generated by the non-vanishing vacuum expectation value of the Higgs 
field (0|$|0) = 771(1)- This mechanism yields on tree-level m w = ^, m u ^ = and 
mu — Af . The renormalised masses on 1-loop level are explicitely given in appendix 
|A] and will be needed later for obtaining the selfconsistent solution. 
In order to be consistent with previous works we choose to scale the quantities in the 
following way: 

mw 



m, 



A a u - m w A" $ - 



Furthermore we use the following representation of the Dirac matrices 

y=(? J I y=( ° ^ -M" 1 ?). (3) 



which allows us to reduce the fermion fields to two components: 

,u,d 3/2 [i>L \ ,u,d 3/2 ( \ 
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With these replacements the new Lagrangian becomes 

v H F ^ Fr + ^p^) f (^) - ^ ( $t$ - 4 ) 2 ) ) > 

with the mass matrix 



(4) 





(5) 

The masses have to be taken from experiment. However, when comparing the masses 
of the different fermions, one finds that the top quark largely dominates the whole 
fermionic mass spectrum. Its energy of ~ 175GeV is much larger than that of any 
other quark, so the other fermion doublets can be considered as massless. In || 
the authors find that the influence of massless fermion doublets is negligibly small 
when compared to a massive doublet, so only the top/bottom doublet needs to be 
considered. However, when applying the usual framework to compute the energy 
spectrum, only doublets that are degenerate in mass can be computed. Consequently 
the top/bottom doublet was approximated by 1.5 fermionic doublets with a degen- 
erate mass of rriT and a corresponding bare value of vp = v u = v&. To estimate the 
effect of nondegenerate fermion masses we developed a perturbative method, which 
can be used to justify this approach. In Appendix || and |L4| we present these results. 
The parameters in this work are chosen to be [15| 

g = 0.67, m w = 83GeV, m T = 175GeV. 

The Higgs mass will be considered as a free parameter. One aim of this work is to 
express the sphaleron energy as a function of the Higgs mass. The relation of the 
bare parameters v F and v H to these physical parameters is discussed below. 
Now we will focus our attention to static solutions of the gauge and Higgs-fields. In 
this case it is useful to use the temporal gauge with A = 0. In this gauge (f|) yields 
the following energy functional: 

£c,ass = ^Jd 3 r Q /-/•;;•' + \ (A*) f (A$) + l/ H ~ 4) 2 ) . (6) 

Configurations with A\ = and $t$ = 4 have -E c iass = which means they are 
vacuum configurations. Due to the transformation properties of (f|) one can trans- 
form this configuration to A\ = iU(x)diW(x) (with U(x) G SU(2)) where the new 
configuration must still have -Eciass = 0, due to gauge invariance. 
The field configurations can furthermore be characterised by the following functional 
(called Chern-Simons-number): 

N ° S= vtrJ + l^aMc) ■ (7) 
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One can show that for configurations of the form Ai = iU(x)diW(x) the Chern- 
Simons-number is an integer value Ncs £ 2- If Acs is not an integer then (||) will yield 
an energy larger than 0. Thus it is impossible to deform vacuum configurations with 
different Ncs continously into each other without violating the boundary condition 

-^class 0- 

Of special interest is the path of lowest energy going from Ncs = to Ncs = 1. It can 
be obtained by minimising the energy functional @ for a given and fixed Acs- This 
path has been constructed in |Tj|. From symmetry arguments the point at Ncs = 1/2 



is the highest point and is called sphaleron. It is a saddle point in the space of Higgs- 
and gauge fields. There is only one negative eigenmode, which is lying in the direction 
of changing Acs- As this point already has a high symmetry in the space of Higgs- 
and gauge-fields it appears logical to choose a spherically symmetric ansatz for the 
spatial dependence of the fields. This ansatz is called hedgehog and has been used 
in a variety of works on this subject [JIB], [TJ, [TjJ]. The Higgs- and gauge-fields are 
represented by 5 functions A(r), B(r), C(r), G(r) and H(r) in the following way: 

1 ~ A (r) sB(r) C(r) 
A"(r) = eaijUj h {b ai - n a ni) h n a rn 



$(r) = 2(H(r)+iG(r)n-r) m. (8) 
In this ansatz the energy (|^) has the following form: 



r 



-) 




2r (H'G - G'H) C - ABGH + 2A (G - H 



+ A 2 + B 2 + ic 2 ) (G 2 + H 2 ) + V \ 2 (G 2 + H 2 - l) . (9) 
For iVcs one obtains the following expression 



2tt Jo 

The fermionic fluctuations have been investigated in || . Here we use only the expres- 
sion for the renormalised fermionic field energy, which is a function of some numerical 
parameters and the physical quantities up, vh and mw- An additional parameter has 
to be introduced, the renormalisation point v Ten . Of course, the final results should 
be independent of z/ ren . Note that for the vp and v H one has to use the renormalised 
quantities on the 1-loop-level (see appendix 0). 
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The numerical parameters, as introduced in are the box size R and the cut-off A. 
Thus the whole fermionic energy at zero temperature becomes 

££ r n m (A, R, v F , u H , v Ten ) = £ sea (. . .) - E div (. . .), (11) 

where the sea energy and the divergent energy are 

^ a (A) = ^r^Trfe^ 2 -e-^ 2 V (12) 



4v^F J\-2 t 3 / 2 
^ div (A) = |5 / (4 n - A 2 ) vl ($t$ _ 4 



rcn 



+^z/^($ t $-4) 2 + 8($ t $-4)JJ. (13) 
The fermionic Hamiltonian is defined to be 

The task of evaluating the energy contribution (|TJ) by diagonalising ([14]) in an ap- 
propriate basis has been described thoroughly in |J and shall not be repeated here. 
Before we further investigate the 1-loop energy expressions we shall examine the 
parameters of the Lagrangian on 1-loop, up, uh and z/ ren . On tree level, the relations 
between these parameters and the physical particle masses are trivial and have been 
given above. But on 1-loop-level these relations become extremely complicated. One 
gets a set of three equations depending on the three parameters. One equation can 
be eliminated to define the mass scale of the system. It appears natural to choose 
the equation for mw- Thus, finally 

m H = rn H (vTen,VF,v H ), 

m F = m F (v iea ,Vp,u H ). (15) 

These system of equations is overdetermined and we have the freedom to choose one 
parameter arbitrarily. Usually one takes 

^ren *^ 

free parameter and examines if the 
final results are independent of z/ ren . The full set of equations ( |15"1) is given in appendix 

E 

Furthermore we need the vacuum expectation value of the Higgs field in the 1-loop 
case. To obtain this value we have to examine the value of the fermionic energy in 
the vacuum which is given by (r is Euler's constant): 

m w N c 



*EL[*3 = 0,* = const] = ^Jd 3 r(-4vL&$ 
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By performing jr^k — one gets the equation for V\\ 



This equation is to be solved numerically for given v F , u H and f ren . It should be 
noted that there is an apparent dependence of Vi on z/ rcn , but by renormalising the 
mass-scale (see above) this dependence is removed. Furthermore it should be noted 
that massless fermions do not influence the value of V\\ since the energy in fll6|) is 
oc up the influence of massless fermions is 0. 

With the help of the expression for the classical energy (the tree level functional) 
(|jD and the 1-loop fluctuations in (|TT|) the correct expression for the 1-loop energy 
functional can be obtained (now we have to use the 1-loop values for up, vu and v\): 

E[A, $] = E c i ass + Ef™ m - (-Edass + -Efcrm) \&$=vf 



9 2 



+E^ m - E- m \^ =vf + ^f(v\ -A) J d?r ($t$ 



After the following replacement 



vi vx 2 2 

2 2 i>i t>i 



one finally gets the following form of (|T§1) 

*[A,»] = I (f^ + CI4 + =^(1-4) M (*** - «)) ■ (19 > 

The task of finding the sphaleron configuration now requires minimising the functional 



19l) at a given value for Nqs- The question is if this value is bound to be Ncs = 1/2 
as it was in the classical case. Clearly the answer is no; since the whole path is not 
symmetric (if the vacuum at iVcs = has an energy of E = then the vacuum at 
JVcs = 1 will have the energy of the created valence fermions which is clearly > 0) 
we would expect the highest point to lie at iV cs > 1/2, but still very close to it. 
In principle, the difference in energy should be less than the energy of the valence 
fermions which is ~ 1% of the sphaleron energy. So it appears to be justified to take 
the minimal value of ([19]) at Nq$ — 1/2 as an approximation for the sphaleron energy 
without introducing larger errors. 
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3 Numerical results 



Now we wish to minimise the functional (|T9|) at the fixed value for Ncs = 1/2. In 
the 1-loop case with the very complicated expression ( |TT|) this is a very complex task 
especially since the condition for Nqs cannot be cast into some explicit form for the 
hedgehog-fields. The problem arises from the fact that the sphaleron configuration is 
rather a saddle point than a minimum. So we have to modify the functional in a way 
that allows us to apply some minimisation procedure. 

The most natural modification involving the boundary condition iVcs = 1/2 is to add 
a term of the form 

a(N cs -l/2) b 

with some parameters a and b. It has turned out that the choices of b = 2 (higher 
values would make the minimum too large) and a = lOOOmw (which is one magnitude 
higher than the sphaleron energy) give good results. So finally one ends up with a 
new functional 

Ex = E X [A, $] + a (N CS [A, $] - 1/2) 2 . (20) 

Ex may either be (|j) or fl!9D since the procedure may also be applied to the classical 
case. 

At this moment we still have a gauge freedom we can use to fix the shape of the 
C-field appearing in the hedgehog ansatz in (H). We choose the following form which 
turned out to be useful for our purposes 

e(«)~»*® 4 (!-i) 4 . 

The quantity s has been set to s = 3m^ which is the size of the other hedgehog 
fields, too. 

In order to do the actual minimisation, the fields A, B, G and H have been discretised. 
For the minimisation procedure, the Powell method (see e.g. jlSfl ) was used and the 
fields are allowed to be modified in the range from r = to r = -R sp hai- -Rsphai is thus 
the maximum size of the sphaleron that can be found with this method. Obviously 
-Rsphai has to be chosen such that 

-Rsphai < R, (21) 

where R is the radius of the box for the Kahana-Ripka basis (see ||). The problem 
is that the symmetries of the Lagrangian (f|) are destroyed when one comes close to 
the point r = R; so a minimisation procedure that searches for a sphaleron in the 
free space (where R — ► oo) will encounter instabilities when it tries to manipulate the 
fields too close to the boundary. This problem is very severe if one uses a minimisation 
procedure which involves computing the functional derivatives of the fields (see |fT9[| ). 
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n a 


n h 


n g 


n h 


Ex/GeV 


31 


21 


21 


35 


12353 


59 


43 


43 


59 


12356 


149 


109 


109 


149 


12356 



Table V. The dependence of the self-consistent solution on the number of lattice 
points. 



mil / GeV 


£ S ,o/GeV 


-E£,mim/GeV 


83 


8804 


8746 


124.5 


9244 


9189 


166 


9497 


9356 


249 


9842 


9839 


415 


10538 


10534 



Table 2: The selfconsistent solution at different Higgs masses. The second column 
shows the initial configuration (the non-selfconsistent Sphaleron) and the third col- 
umn the final configuration after the self-consistent minimisation. One can see that 
the differences are minimal and the total energy is well approximated by the initial 
configuration. 

However, since our method does not involve derivatives it should be stable and reliable 
if we ensure that ( ffi] ) is satisfied. 

Before performing this examination the lattice spacing should be examined first; the 
number of lattice points of the different fields are presented in table |l[ The bare 
parameters of the Lagrangian have been chosen to be vu = 0.8, vp = 2.1 and u ren = 
2.0. We find that the first row in the table is fine enough to represent the fields. 
In order to determine an appropriate value for i? sp hai, we examined the energy func- 
tional ([19]). We found that the total energy does not depend on the choice of -R sp h a i/-R, 
and it is convenient to use -R sp hai = 0.70i? as a good compromise between accuracy 
and computational effort of the problem. 

Now we can examine three different solutions: the tree-level sphaleron (the sphaleron 
solution of (0) considering only the classical energy of the boson field), the one-loop 
sphaleron (the field configuration of the tree- level sphaleron evaluated with (|i~9|) ) and 
the selfconsistent sphaleron (the sphaleron solution of ([19])). 

By computing the energy of these solutions and comparing the deviations one can 
make a statement about the significance of self-consistency in the one-loop approxi- 
mation of the sphaleron. The actual results for different values of the Higgs mass are 
listed in table ^. 
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One should keep in mind that the errors of the selfconsistent calculation are about 
1% due to the inaccuracy of the minimisation process. However the influence of the 
self-consistency is negligible, and the simplest way to compute the one-loop sphaleron 
energy seems to be to minimise the classical part of ([19]) and then add the fermionic 
contribution. Of course one has to use the one-loop renormalised parameters in 
evaluating fllTf ) and in computing the final energy. 

Sphaleron energy 



1Q5QQ 
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Figure 1: Comparison of the sphaleron energies for different values of the Higgs mass. 
The dots represent the energies of the classical solution, the crosses represent the 
energies of the one-loop solutions. The circles represent the results based on the 
calculations in || using the renormalization from [|TIJ . 



With this recipe at hand we can compare the one-loop and the tree-level energies of 
the sphaleron for a wider range of Higgs masses. This has been done in figure |l|. One 
should remember that the vacuum configuration is no longer stable when the Higgs 
mass becomes too small; this phenomenon has also been used to find a lower bound 
for the Higgs mass (see e.g. |20|). However, in this work it only causes that equation 
(|17D cannot be solved numerically and we cannot compute the one-loop solution. 
It is clearly visible that the shape of the configuration is maintained in comparison 
to the classical curve. But the overall energy has increased by an amount of about 
2% — 3%. From figure [I] we can thus conclude that the influence of the one-loop 
fluctuations is only a few percent in comparison to the tree level result. 
When comparing these results with the results from previous computations || |HJ , we 
find that the previous results deviate significantly from the current results. The old 
curve (represented by circles in [l]) is much too large for small values of the Higgs mass 
and even exhibits a local minimum at about m# ~ m\y- This is a strong indicator 
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for the incomplete nature of the renormalization scheme used. 



4 Conclusion 

In the last section we have found that the correction to the energy of the sphaleron 
configuration due to fermionic contributions is very small, only about 2% — 3%. Fur- 
thermore in computing this effect the influence of self-consistency turned out to be 
negligibly small. Since the fermionic contribution is only small in principle the shape 
of the dependency of the sphaleron energy on the Higgs mass did not change - it is 
only shifted by a small amount upwards. This puts the results in || into doubt, since 
they found due to the incomplete renormalization a strong increase in the sphaleron 
energy even at zero temperature. 

Still, an interesting question is whether the results of this paper can be generalised to 
the case of non-zero temperature. In this case the correction generated by the fermion 
determinant is not much smaller than the classical sphaleron energy. Although it 
appears possible that now self-consistency becomes important, it has been shown in 
H that the fermionic contribution can be split up in two parts: the major part has 
the same shape as the classical energy and the remaining part can be treated as a 
small correction. Now one can include the first part of the fermionic energy into 
the classical one and compute the sphaleron energy again (practically this means one 
has to adjust the bare-parameters again and then proceed in a similar way as in this 
paper). But it was shown in that then the remaining fermionic contribution again 
has the effect of a small correction to the overall energy. So in this case the energy of 
the whole configuration should not be significantly influenced by self-consistency. So 
we argue that also in the case of finite temperature it is unlikely that self-consistency 
introduces any changes in the energy. 



However, the results obtained in |I0| appear to suffer from the same shortcomings 
(incomplete renormalization) as those in || when it comes to the absolute value of 
the energy. Thus the conclusions drawn for the upper limit of the Higgs mass seem 
to require a correction that would lower the value and thus move it closer to the one 
cited by Shaposhnikov in JT2 . 



Acknowledgement: The authors wish to thank D. Diakonov, P. Pobylitsa, M. Polya- 
kov, J. Schaldach, P. Sieber and C. Weiss for their valuable input and discussion. 
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A Mass renormalisation 



Here we want to present the shape of the equations fll5l) . These were calculated in 
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On tree level we have 

mw = 83GeV, 
m H = v H m w , 

m F = v F mw (22) 

In the case of one-loop fluctuations we find the following set of equations to determine 
the bare parameters: 



= -rn 2 H + ^{3vl 



rcn , 



rcn 



2 



.2 , V l 



-m w + 4 

92 ( Nc ir? F> | -n , 



\ \ ' run ^ ' ftu I / \ \ 1(, rmi ' ' 



ren ren/ \ \ rcn " ren, 

2 / 2 22\ / 2 2 2 ^ 



^ren V ^ron ' 4l/ ren / \ 4z/ren2 4z/ rc 

x -C-G 



z/ 2 '4z/ 2 

i ren ren , 



X ( iV + ^ - ^ (-C - G (^,0 



a 1 9 / l "ren O I I ,,2 ' 

Z 7 \ Z \ \ 'rcn / / 

= \v>l-™ F , (23) 
where we have used the following definitions: 

-1/2 



G (x, y) = J da In (^x (^a 2 — + y — ie^j 

r 1 / 2 ( ( 1 

G' (x, y) = J da a 2 In yx ya 2 — — 



is 
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B Nondegenerate Fermion Masses 



As mentioned in the main text, it is not possible in our approach to treat nondegen- 
erate fermion masses exactly. So we introduced a common fermion mass v F and a 
degenerate mass matrix 



Mn 




(24) 



to compute the results presented here. In this appendix we try to estimate the error 
of that approximation by treating nondegenerate fermion masses perturbatively. In 
order to restore the original mass matrix M, (|5]), one can write 



Av 

M = M 1 + r 3 

Up 



and Av 



v„ 



(25) 



One can see 



where we introduced the parameters Vp ^ ^ 

that the matrix r 3 spoils the spherical symmetry of our fermionic Hamiltonian (13) 
so that a direct diagonalisation can not be performed. But for small mass differences 
one can now expand M in a power series with respect to Avj v F . Of course, in nature 
this parameter for the massive fermion doublet is of order 0(1), so it is clear, that 
it is not possible to compute good estimates with this ansatz. On the other hand, if 
we are only interested in the qualitative behaviour of the fermionic energy, it can be 
used to see how the energy changes if one allows different fermion masses within one 
doublet. The fermionic 1-loop energy can then be evaluated as a Taylor series in the 
parameter Av: 



EFerm(Av) = Ep erm (0) + 



83 



Ferm 



dAv 



■ Av + 



1 d 2 Ep erm 



2 dAv 2 



Az/ 2 + 0(Az/ 3 ).(26) 



In this series the linear term vanishes, because in a pure SU(2) gauge theory up- and 
down-components in one doublet cannot be distinguished if they are degenerate in 
mass. Therefore we compute numerically the first and third term of this series. To 
compare our result to the approximation we proceed in the following way: We start 
with degenerate fermion masses with v F = v t and count only n = 3/2 massive fermion 
doublets, so that in the end three top quarks masses are counted. This is exactly the 
approximation which was used before. Then we introduce small mass differences Av, 
but keep the value of v t fixed. To ensure that we have the correct "mass content" we 
count 



n 



2 v t - Av 



(27) 



massive fermion doublets, in order to preserve the sum of three top quark masses in 
the fermionic energy. 
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By increasing Au we eventually reach the "physical" point with Au = u t /2 and n = 3. 
This is plotted in figure |||. This method will work well in the region of small Au, 
since perturbation theory is well justified here. We can clearly see, how the fermionic 
energy drops if one introduces mass splittings. From this point of view we can at 
least say that our approximation is a good upper limit for the energy arising from 
fermionic fluctuations. 



5.0 



EFerm/mW 




2.0- 



1.0- 

0.0 I 1 1 1 1 1 1 1 1 1 

0.0 0.1 0.2 0.3 0.4 0.5 

Figure 2: In this picture one can see the fermionic energy as a function of the mass 
splitting Au. The dashed line is the first term of the series fl26|) and the dotted line 
shows the third term. The solid line is the sum of the two contributions and therefore 
the whole fermionic energy up to Az/ 4 , since all odd powers vanish. One can see that 
even a small mass difference lowers the fermionic energy significantly. 
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